Energy dependent counting statistics in diffusive superconducting tunnel junctions. 
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We present an investigation of the energy dependence of the full charge counting statistics in 
diffusive normal-insulating-normal-insulating-superconducting junctions. It is found that the current 
in general is transported via a correlated transfer of pairs of electrons. Only in the case of strongly 
asymmetric tunnel barriers or energies much larger than the Thouless energy is the pair transfer 
uncorrelated. The second cumulant, the noise, is found to depend strongly on the applied voltage 
and temperature. For a junction resistance dominated by the tunnel barrier to the normal reservoir, 
the differential shot noise shows a double peak feature at voltages of the order of the Thouless 
energy, a signature of an ensemble averaged electron-hole resonance. 
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The mechanism for current transport across a nor- 
mal conductor-superconducting (NS) interface, at ener- 
gies below the superconducting gap, is Andreev reflec- 
tion. An electron incident from the normal conductor 
is retroreflected into a hole, and a charge of 2e is trans- 
ported into the superconductor. The fact that electrons 
are transfered across the NS-interface in pairs has re-, 
cently attracted a lot of interest to noise in NS-systems,EJ 
since the noise can provide information about the charge 
transfer mechanism. 

In some NS-systems, the shot noise is doubled com- 
pared to the noise in the corresponding normal system. 
In a normal-insulating-superconducting (NIS) junction, 
the pairs tunnel with a small probability between the 
normal and the superconducting reservoirs. This gives 
rise to a Poisson noise of pairs of electrons™|«dhch is thus 
twice as large as standard Poisson noise.Hau Such pair 
Poissonian noise was recently observed experimentally!! 

Also in diffusive wire NS-junctions, a douMLag of 
the shot noise has beenjpredicted theoreticallyB'Qu and 
observed experimentallyQtrt3 for applied voltages much 
smaller as well as much larger than the Thouless energy. 
However, for intermediate voltages, of the order of the 
Thouless energy, .the shot noise deviates from the dou- 
ble normal value,l3Eil a behavior related to the induced 
proximity effect in the normal conductor. 

In general, in mesoscopic NS-systems, there is no 
simple relation between the shot noise in the normal 
and the superconducting state. This has been demon- 
strated for a wide variety n aL multi-mode systems such 
as diffusive NIS-iunctions,crE3 chaotic-dot superconduc- 
tor junctions^ 3 b 14 1, multiterminal superconducting .tup-, 
nel junctions^ said wide ballistic double-barrierE3t£l 
and normal slabO superconductor junctions, as well 
as few- m ode . supe rconducting systems containing beam- 
splitters,Hli§E30 disordered tunnel junctions^ and dou- 
ble barrier junctions.E3 

To obtain complete knowledge of the charge trans- 
fer in NS-systems, it is thus not sufficient to study 
the shot noise. Instead, one has to investigate the full 
charge counting statistics.Ej This was originally done by 
Muzykanskii and Khmelnitskiip who showed that the 



statistics of the charge transfer across a NS-interface with 
arbitrary transparency, can in general be described as 
a correlated transfer of pairs of electrons. Only in the 
limit of a low transparency NIS junction, the transfer 
of pairs is uncorrelated, Poissonk^pJJaing a recently de- 
veloped circuit theory approachLHl3o'E£l to full counting 
statistics, the charge transfer statistics has been stud- 
ied in several systems containing superconductors, such 
as short SNS-j unctions £3 diffusive NS-j unctions ,0 diffu- 
sive superconducting tunnel junctionals and chaotic dot- 
superconducing junctionsEH However, in these works, the 
main focus has been on the low voltage- and temperature 
properties, the full energy dependence of the counting 
statistics was only analyzed numerically in Ref . j?J . 

In this paper we investigate the energy dependent 
full counting statistics of a diffusive normal-insulating- 
normal-insulating-superconducting {NI\N' I2S) junc- 
tion. This system is interestiutiot-several reasons. Using 
the circuit theory approach JjtLaEa an analytical expres- 
sion for the full voltage and temperature dependence of 
the counting statistics can be derived. This allows us to 
investigate in detail the energy dependence of the charge 
transfer mechanism. We focus on the role of the proxim- 
ity effect, which is of particular interest, since it gives rise 
to a large, energy dependent modification of the trans- 
port properties. Moreover, the conductance inJSI I\N' I2S 
systems was recently studied experimentallyJ23 and very 
good agreement was found between the experiment and 
the quasiclassical Greens function theoryO on which the 
circuit theory is based. This makes a detailed analysis of 
the experimentally accessible noise of interest. 

We find that the current is in general transported via 
a correlated transfer of pairs of electrons. However, in 
the cases where the resistance is dominated by the tun- 
nel barrier 12, at the N'^S interface, or for energies well 
above the Thouless energy, the pair transfer is uncorre- 
lated. These cases correspond to the limit of weak prox- 
imity effect in the normal region N'. In the opposite 
regime of strong proximity effect in AT', when the resis- 
tance is dominated by the tunnel barrier I\ at the NI2N' 
interface, a proximity gap opens up in the ./V'-region. For 
energies below the induced proximity gap, the junction is 



2 



effectively a NIS-junction and the pair transfer is again 
uncorrelated. 

The second cumulant, the noise, is found to depend 
strongly on the applied voltage and temperature. For 
a junction resistance dominated by the tunnel barrier 
I\, at the NI\N' interface, and a temperature well be- 
low the Thouless energy, the differential noise shows a 
double-peak feature at voltages of the order of the Thou- 
less energy. We show that this double-peak behavior can 
be attributed to an ensemble averaged electron-hole res- 
onance. 

The paper is organized as follows. We first present 
the circuit theory model and the derivation of the full 
charge counting statistics. Thereafter, the dependence of 
the counting statistics on energy and tunnel barrier (7i 
and I2) conductances is analyzed. We then investigate 
the temperature and voltage dependence of the noise in 
detail. 



I. MODEL AND THEORY 

We consider a diffusive normal-insulating-normal- 
insulating-superconducting (NIiN'I 2 S) junction (see 
Fig. |l|), where both the resistances of the insulating 
tunnel barriers (ii and I2) are much larger than the 
resistance of the normal diffusive region (N 1 ) between 
the barriers. The short diffusive region N', which does 
not contribute to the junction resistance, leads to com- 
plete isotropization of the system parameters between 
the barrierSjtJ This makes the junction effectively zero 
dimensionality and the exact shape and size of the diffu- 
sive region is thus not important. 
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iVi(2) is the number of transport modes at the barrier 
1(2) and r 1 ( 2 ), n is the transparency of mode n. We con- 
sider the limit of large number of modes, ^Vi(a) ^> 1 an d 
small transparency of each mode, r 1 ( 2 ),n ^ 1- How- 
ever, the conductances g\ and g 2 are much larger than 

the conductance quanta 2e 2 /h, i.e. 53n=? ) ^1(2), n ^ 1- 
Under these conditions we can neglect Coulomb blockade 
effects as well as weak localization like corrections to the 
cumulants of the current. 

It is further assumed that the dwell time of the par- 
ticles in the the normal region N is much shorter than 
the inelastic scattering- and phase breaking times. We 
are also only interested in zero frequency properties.t3 
With these assumptions, we can analyze the full count- 
ing statistics of the junction within the framework of the 
circuit theory, recently, developed by NazarovElEij and 
Belzig and NazarovflE3 Very recently, a detailed discus- 
sion of the circuit theory was given.Ea We note that the 
zero energy counting statistics was studied for the same 
junction in Ref. (lq| , here we study the full energy de- 
pendence. 

In the circuit theory, the mesoscopic conductor is de- 
scribed in terms of nodes connected via various circuit 
elements, just as in an ordinary electric circuit theory. 
Due to the effective zero-dimensionality of the junction, 
the circuit theory representation (see Fig. [l]) only con- 
tains three nodes, the two reservoirs (N and S) and the 
isotropic normal region (N') between the tunnel barriers, 
and two circuit elements, the tunnel barriers (I\ and I 2 ). 

However, instead of characterizing each node with an 
electric potential, as is done in an ordinary electrical cir- 
cuit theory, they are characterized with a 4 x 4 matrix 
Green's function (in Keldysh-Nambu space). Here we 
have the three node Green's functions Gn, Gs and G, as 
shown in Fig. pi They obey the normalization condition 
G 2 N = G| = Gr = I. Throughout the paper we use the 
convention with "check" , 7 denoting 4x4 matrices and 
"hat" ? denoting 2x2 matrices. 

The matrices of the reservoir nodes are known. In the 
normal reservoir we have 
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FIG. 1: The JV/iiV'^S'-junction. Upper figure: A normal 
(N) and a superconducting (S) reservoir are connected via 
a diffusive normal metal (N' , shaded) contact. At the NN' 
and N'S interface, there are tunnel barriers I\ and I2 with 
conductance g\ and Q2 respectively. Lower figure: The circuit 
theory representation of the junction, with the matrix Greens 
functions Gn, Gs and G and the counting fields xn and xs 
shown. 

The conductances gi and 172 of the two tunnel barriers 
can be written as gi( 2 ) — (2e 2 //i) J2n=i } ^i(2),n, where 
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where f±(E) = f(E±eV), with f(E) = (exp[.E/fcT] + 
1) _1 the Fermi distribution function. Here \n is a count- 
ing field, "counting" the electrons passing in and out of 
the normal reservoir and a z is a Pauli matrix. In the 
superconducting reservoir we have 



G° s 



K, 



Rs K s 
A s 



R(A) t 



_ ( 9R(A) fR(A) 



-9r(a) 



s — \ f 

JR(A) 

(Rs-A s )[f(E)-f(-E)}, 



(2) 



where A is the superconducting gap and /r(a) — «A[(_B± 
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ie) 2 -A 2 } 1 / 2 a,ndg R ( A ) = [l-/^^)] 1 ^ 2 are Green's func- 
tions characterizing the superconducting reservoir (e is 
an infinitesimal positive number). The counting field xs 
of the superconducting reservoir is introduced for nota- 
tional clarity, due to current conservation it is in principle 
not needed. 

The matrix G of the node between the tunnel barriers 
is to be determined. It is found from a conservation law 
for matrix currents (a matrix Kirchoff's rule), as 



h + h + Ie = 
where the matrix currents are given by 

h 



(3) 



Ie — j- [Ge , G] , Ge = ( ^ 



(4) 



where [A, B] = AB — BA and 8 is the mean level spac- 
ing in the A^'-region divided by the conductance quanta 
2e 2 /h. The Thouless energy is given by Em = (gi + 
92)5/4:, i.e. it depends only on the escape time through 
the tunnel barriers, not on the diffusion time through 
the iV'-region. The term Ie leads to a modification of 
the proximity effect at finite energies, and eventually a 
suppression for E 3> E^k- 

From Eq. (||) and the condition G 2 = 1 , we can deter- 
mine the missing node matrix G(xn, Xs) as a function of 
the counting fields. Following the lines of Ref. |l5| , not- 
ing that Eq. (||) can be written in the form [G, G x ] = 0, 
we find the physically relevant solutions 



G = G X (G 2 X )- 1/2 
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Knowing G(xn,Xs), we can find the cumulant gener- 
ating function F (xn,Xs)- It is given from the general 
relation between the matrix current and the generating 
functionH3 
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where to is the measurement time. Using that under the 
trace, 

djGl) 1 / 2 dG x v ,^2\-i/2 i \~ ifi 
-5 - ^ G X (G X ) — -gti 2 )[TK, G N / S )\G 

= -^9i(2)Tk[G N ( S )-, G] = -f K h(2), (7) 

we get the counting field dependent part of the cumulant 
generating function 
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We note that the standard definition of a generating func- 
tion demands F(xn — 0, xs — 0) = 0. This implies 
that a counting field independent constant — tr[G 2 (xN = 
0;XS = 0)] 1 / 2 should in principle be added to the inte- 
grand in Eq. (||). However, such a constant does not alter 
the cumulants, and below we will for simplicity omit it. 

The probability P(Q) that Q electron charges have 
been transported through the junction in time to is given 
by 



P(Q) = 



dx 



N 



2n 



-iXnQ-F{xN ,xs=0) 



(9) 



Eqs. (||) - (^|) completely determine the full charge count- 
ing statistics for arbitrary temperatures and voltages. 
For our purposes, it is however more convenient to focus 
on the properties of the cumulant generating function in 
Eq. (|^). From the cumulant generating function we can 
derive all cumulants of the current by repeatedly taking 
the derivative with respect to the counting field xn (or 
Xs)- For the two first cumulants, current and noise, we 
have (evaluated at the iV-contact), 
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Xn=Xs=0 

and similarly for higher order cumulants. 

II. FULL COUNTING STATISTICS. 

In what follows we focus on the limit kT, eV, Etk <C A. 
In this limit, the physics is dominated by the induced 
proximity effect in the diffusive region N' between the 
two barriers .L3 The proximity effect affects the transport 
properties on the energy scale Etk- It can on a micro- 
scopic level be described as resulting from coherence be- 
tween the two quasiparticles in the Andreev reflection, 
e.g. an incident electron and a retrorefiected hole. An in- 
duced proximity effect leads in general to large modifica- 
tions of the transport properties, i.e. modifications scal- 
ing with the system size. We note that this is in strong 
contrast to the situation in normal systems, where phase 
coherence only gives rise to small quantum corrections. □ 

However, being a coherent phenomenon, the proximity 
effect is also sensitive to modifications of the electron and 
hole phases. At finite energies, the electron and hole pick 
up different phases and at E ^> Eth, the electron-hole 
coherence is destroyed and consequently, the proximity 
effect is suppressed. Moreover, a magnetic flux in the N'- 
rcgion large enough to break time reversal symmetry, also 
suppresses the proximity effect, since the electron and the 
hole pick up different phases in a finite magnetic field. 
Here we however only consider the case with negligibly 
small magnetic fields. 



4 



Formally, at energies well below the superconducting 
gap A, the matrix Greens function G s in Eq. (0) simpli- 
fies considerably, since R$ = Ag — a x and consequently 
Ks = 0. We then find F(xn,Xs) in Eq- (H) by diago- 
nalizing (G 2 ) 1 / 2 , giving 



F(xn,Xs) 
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where we have introduced 



2 1 2 
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A = (l- e 2i (x«-«)^ 

+ (l-e- 2 ^"-«))/_(l-/+). (13) 

We note that at E = 0, we recover the generating func- 
tion of Rcf. |fl5|| . The terms exp(±2z[xjv — Xs]) show 
that electrons are transported in pairs out of (+) and 
into (— ) the superconductor and into (+) and out of (— ) 
the normal reservoir. This pair transport follows from 
the fact that the charge transfer mechanism across the 
iVS'-interface is Andreev reflection. The form of the cu- 
mulant generating function, with the double square roots, 
tells that the pair transfer is highly correlated. There 
are however several limits where the generating function 
is simply proportional to A, describing an uncorrelated, 
(generalized) Poissonian, transfer of pairs of electrons (as 
follows from Eq. |9|). 

In the limit of high energy, E 3> Eth , for arbitrary 31 
and g 2 , the integrand of the generating function in Eq. 
( |l2| ) takes on the form 



9F(xn,Xs) 
dE 



g\g\ 



4e 2 (31 + <? 2 ) 2 
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For these energies, the probability of pair tunneling is 
strongly suppressed and the pairs are thus emitted in a 
Poisson process. This is a consequence of the different 
phases picked up by the electrons and holes, ~ zLE/Exh, 
at finite energies. At E ^> Eth the electrons and holes 
loose their coherence and the pair tunneling probability 
is suppressed. 

In the limit of a dominating coupling of the diffusive 
region N' to the normal reservoir, 31 ^> g 2 , the integrand 
of the generating function becomes 



dF(xN,Xs) h 



9l 



1 



dE 



4e 2 gi 1 + {E/E Th ) 2 



A. 



(15) 



Thus, in this limit, the pair transfer is also uncorre- 
lated. In this case, the barrier I\ has a negligible re- 
sistance compared to I 2 and the junction is effectively a 
A^S-j unction. As discussed above, in this case the pair 



transfer is uncorrelated. We note that for high energies 
E > E Th , it coincides with Eq. (|l4|), (for 31 > g 2 ). 

We also noteta that for the limit 31 <C 32, strong cou- 
pling of N' to the superconductor and for the additional 
limit of low energies, E <C Exh, the pair transfer is un- 
correlated. In this case, the integrand of the generating 
function is simply 



dF( X N,Xs) 
dE 



4e 2 



.92 



A 



(16) 



The reason for the uncorrelated pair transfer in this case 
is the following: In the limit 31/32 — > 0, a proximity gap 
opens up in N' and the junction becomes an effective 
NI\S' junction (5" denoting the gapped normal region 
between the barriers). For energies below the induced 
gap, the pair transfer is thus Poissonian, as in a standard 
AT75-junctions.M 

Making a connection to the proximity effect, it can be 
noted that the in two cases, Eq. ( |i"4| ) and Eq. (|i"5|), the 
strongly suppressed proximity effect in N 1 is responsible 
for the uncorrelated pair transport. However, in the last 
case, Eq. (|l6|), the Poissonian pair transfer is due to the 
induced proximity gap in N' . 

It is instructive to also study qualitatively the full 
counting statistics, starting from the scattering point of 
view in Ref . [^) . The cumulant generating function, gen- 
eralized to many modes, is 

F(xn, Xs) = jJ dE M 1 + Reh(E)A)) (17) 

where the R e h(E)'s are the Andreev reflection eigenval- 
ues, i.e. the real and energy dependent eigenvalues of the 
hermitian electron- hole scattering matrix product S eh S e h 
(see e.g. Ref. H). Here (..) denotes ensemble average 
over impurity configurations and A is the same as in Eq. 
(|l3|). For a general function h(R e h), the ensemble average 
can be written as an integration J dR e hp(R e h)h(R e h), 
where p(R e h) is the distribution of Andreev reflection 
eigenvalues (the corresponding quantity in NS-systems 
to the transmission eigenvalue distribution in normal sys- 
tems). 

Expanding the logarithm in Eq. ( |l7| ) , makes it is clear 
that the limits where the charge transfer is uncorrelated, 
Eqs. ( |l4"| ) to ([l6|), corresponds to the limit where (R e h) ^> 
(R 2 h ) (and consequently (R e h) 3> (R"h) f° r au high 01 ' 
moments n > 2 since < R e h < 1). 

From this we can draw the following conclusions. The 
limits of uncorrelated pair transfer corresponds to the 
cases where there are no or very few "open" Andreev 
channels, i.e with Andreev reflection probability close to 
unity (leading to (R e h) ^> {R 2 e h))- ^ n the case where there 
are more open channels, the various moments (R" h ) are 
of the same order of magnitude, resulting e.g. in a noise 
below the Poissonian noise of uncorrelated pair transfer 
across the NS-interface, further discussed below. As is 
also further discussed below, this Andreev channel pic- 
ture provides a simple explanation for the voltage and 
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temperature dependence of the noise as well as the cur- 
rent. 



III. NOISE. 

The rest of the paper is devoted to a detailed anal- 
ysis of the properties of the the second cumulant, the 
noise, which is within reach of existing experimental tech- 
niques. We note that the cumulant generating function 
F(xNi Xs) m Eq. ( p^ ) can be written as g\ (or 172) times 
a function which only contains the ratio (71/32, he. all 
transport properties studied below, (normalized) depend 
only on the ratio of the conductances. 

For the benefit of the following discussion, we first 
study the current, given from Eq. @ and ©, 
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The current has. been studied both theoreticall and 
experimentally!] here we briefly summarize the findings. 
The differential conductance G = dl/dV at zero temper- 
ature is plotted as a function of eV/Exh in Fig. |^, for 
different ratios 171/92 ■ The conductance for eV <C Eth is 
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FIG. 2: The differential conductance dl/dV, normalized with 
the zero voltage conductance Gns, as a function of eV/ Eth- 
The values for the ratios 31/(72 = 5 and the corresponding 
line types are shown in the legend. The temperature kT <C 
eV. For g? — 4gi, the conductance shows a strong peak at 
eV ~ Erh, a signature of an ensemble averaged electron-hole 
resonance. 



given by 



G(eV < E Th ) = G 



NS 



1 1 

9x92 



{9l + 9lf 2 



(19) 



In the opposite limit, eV 3> E^u, the conduc- 
tance decays as a power law with applied voltage, as 
(E th /eV) 2 g 2 g 1 /(g 1 +g 2 ) 2 . 



The behavior of the conductance for voltages of the 
order of Erh depends strongly on the ratio 31/32- For 
9i 3* 92, i-e. for a strong coupling of the normal region 
N' to the normal reservoir, the conductance decreases 
monotonically with voltage as (gf/5i)/[l + (eV/ETh) 2 }- 
In the opposite regime, 32 2> 31, the normal region is 
strongly coupled to the superconducting reservoir. The 
conductance has a strong peak at voltages eV ~ Eth> 
which can be seen as an ensemble averaged electron-hole 
resonance. In other words, for energies E ~ Exh the 
distribution of Andreev reflection eigenvalues, p(R e h), is 
shifted towards more "open" channels. 

As pointed out above, in the limit of .91/32 ™> 0, i.e. 
when decoupling the normal reservoir, a proximity in- 
duced gap Eth opens up in the spectrum in the nor- 
mal region N'. In this limit, the conductance shows a 
singularity at eV = Exh, similar to the standard NIS- 
tunneling conductance, but with the induced gap Eth 
instead of the superconducting gap A. 




kT/Exh 

FIG. 3: The zero voltage conductance normalized with Gns, 
as a function of temperature. The values for (71/02 and the 
corresponding line-types are the same as in Fig. The con- 
ductance shows a temperature dependence which is qualita- 
tively similar to the voltage dependence of the conductance 
in Fig. |i| The zero voltage conductance is related to the 
equilibrium noise as S eq {T) = 4ekTG(eV = 0,T). 



The temperature dependence of the zero voltage con- 
ductance, G(T), is shown in Fig. ^. The temperature 
dependence is qualitatively similar to the voltage depen- 
dence, shown in Fig. pi d ue to the form of the spectral 
current density in Eq. (|18[) . We note that for kT <C Erh 
we have G — Gns from Eq. (^9|). In the opposite limit, 
kT ^> ETh the conductance decays as 1/kT. In the inter- 
mediate regime, the behavior of the conductance depends 
strongly on the ratio gi/g 2 . For 31 3> (72, it decays mono- 
tonically with temperature, while in the opposite limit, 
92 ^ Si) it has a peak at temperatures kT of order Exh- 

We then turn to the noise, found from Eq. ( pT| ) and 
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(H) to be 
S = J dE 



V« 2 + P 2 (a + V« 2 + P 2 ) 
[/+(! - /-) + /-(I - /+)] (a + 

-tf- - /+)2 ^rt ( 2a + 3 v / ^ J +^)| ■ (20) 



First we consider the equilibrium noise, for eV <C 
ETh, kT. In this limit we have /_ — /+ = and the 
noise is 
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dE 



4^2/(1 - j)g\gl 



y/a 2 + /3 2 \Ja + y/a 2 +7^ 



(21) 



Noting that we have 2/(1 - /) = kT[d(J- - 
f+)/dV} eV= a = -2kT(df/dE) we see from Eq. @ and 
@ that 



S eq (T) = AkTG(eV = 0,T), 



(22) 



for arbitrary temperature T. This is just what we expect 
from the fluctuation-dissipation theorem. The temper- 
ature dependence of the noise, normalized with 4kT, is 
thus just the temperature dependence of the zero voltage 
conductance, plotted in Fig. (||). 

At kT < E Th we find just S eq (kT < E Th ) = 
AkTG^s- ln the opposite limit, kT 3> Eth the equi- 
librium noise tends towards a constant, temperature in- 
dependent value 



S eq {kT^>E> 



Th) 
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Th 



(.9i + g2)vgj+gl 



(23) 



In the intermediate temperature regime, the noise, as 
the conductance in Fig. 0, depends strongly on the ratio 
9\l 92- 

The voltage dependence of the noise is qualitatively 
different from the temperature dependence. In the limit 
kT <C eV, Eth, the shot noise regime, we have 
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Here we focus on the differential noise, dS/dV, which is 
plotted as a function of voltage for different ratios 31/32 
in Fig. §. 



In the low vo' 
noise is given byt 



limit, eV <C Exh, the differential 
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The Fano factor, P/vs/2eGjvs, in this limit, is thus 
(2 — 5g 2 g 2 /[gl + g 2 ] 2 ). It varies from 2 in the limits 




FIG. 4: The differential shot noise dS/dV as a function of 
applied voltage, normalized with the zero voltage value. The 
values for <?i lg% and the corresponding line- types are the same 
as in Fig. 0. The temperature kT <C E T h- The double- 
peak behavior at eV ~ Em, for <?2 = igi, is a signature of 
the electron-hole resonance, responsible for the peak in the 
conductance at eV ~ Eth in Fig. [| 

3i 52 and 32 ^ 3i, to 3/4 forEl 31 — g 2 - A Fano fac- 
tor 2 is an indication of uncorrelated emission of pairs of 
electrons into or out of the superconductor. As pointed 
out in connection to Eq. ( |l5| ) and (|l6|), in these limits, 
for eV <C E T h and g\ » 32 or .92 > gi, the generat, 
ing function F(xn,Xs) m Eq- (0) is proportional tcO 
exp(2i[xN — Xs])i showing that the charge transfer is in- 
deed an uncorrelated emission (for negative voltages) of 
pairs of electrons from the superconductor. 

In the high voltage limit, eV 3> Eth, the differential 
noise decreases with voltage as 



edV \eV 



9%9i 



{91+92) 



(26) 



i.e. the same power law behavior as the conductance, 
which is clear directly from the generating function in 
Eq. ©■ 

At intermediate voltages, the noise depends strongly 
on the ratio 31/32- For a dominating coupling of the nor- 
mal region to the normal reservoir, 31 ^> 32, the differen- 
tial noise decreases monotonically with applied voltage, 
as 



1 dS _ gj 1 

edV 3! 1 + (eV/E Th ) 2 ' 



(27) 



This can also be obtained directly from the generat- 
ing function in Eq. (|l5|). We note that this behavior 
is very sisailar to what was found for a diffusive NIS- 
junctionjlj a diffusive normal region of non-negligible re- 
sistance, connected to a superconducting reservoir via a 
tunnel barrier. 

In the opposite limit, gi 3> 32, the noise shows a dou- 
ble peak structure around eV ~ Exh, the voltage where 
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the conductance shows a single peak (see Fig. ||). Such 
a double peak behavior has been discussed before iiMm- 
gle mode junctions with sharp Andreev resonancesJla In 
single mode junctions, the double peak behavior follows 
directly from the fact that if the Andreev reflection prob- 
ability R e h has a resonance at some energy Eq, the noise, 
~ -R e fr(l — Reh), shows a double peak around Eq. In- 
terestingly, as is clear from Fig. [|, in the NIiN'I^S, 
this double-peak behavior survives the ensemble average. 
This indicates that the whole distribution p(R e h) of An- 
dreev reflection eigenvalues is shifted from being domi- 
nated by "closed" towards being dominated by "open" 
channels at resonance E ~ Eth- 

We note that such a double peak behavior can also be 
found in normal single-mode double-barrier junctions, it 
is not an effect of the superconductivity (see Ref. [jy). 
However, in normal, many-mode systems, such a be- 
haviour is, to leading order in number of modes, washed 
out when performing an ensemble average. 



> 

> 

5j 1 

T3 



V ' i 

\ / \ / / 
A > ,' / 



1 2 

eV/E Th 

FIG. 5: The differential Fano factor [dS(V) /dI(V)]/2e as a 
function of applied voltage. The values for gi/02 and the 
corresponding line-types are the same as in Fig. 0. The tem- 
perature kT <C E T h- The dip at eV ~ E T h, for 32 = 4<j>i, is 
a signature of the electron-hole resonance, responsible for the 
peak in the conductance at eV ~ Eth in Fig. |j| 

Additional insight can be obtained by studying the dif- 
ferential Fano factor, [dS(V)/dI(V)}/2e. In the differen- 
tial Fano factor, the voltage dependence of the conduc- 
tance is canceled out from the noise, and we obtain the 
"bare" voltage-dependence of the noise. We have 



1 dS(V) 
2e dI(V) 



3^ 



/3 2 ) 



(a 2 + /3 2 )( 



a + \fo? 



(28) 



This quantity is plotted in Fig. g. We note that 
the differential Fano factor shows a much stronger de- 
pendence—^ .the., applied voltage, compared to earlier 
studiedDBE3e30 diffusive NS and NIS junctions. The 
resonant double-peak in the noise, for gi 3> g\ is not 
visible in the differential Fano factor, instead it shows 
a dip at eV ~ -Erfc- This is in accordance with our 
explanation of the double peak-behavior of the noise 
above, since, again making the comparison to a single 
mode junction, the differential Fano factor is the differ- 
ential shot noise divided by the differential conductance, 
— [R e h(l - Reh)]/R e h = 1 - Reh, which thus has a dip 
at the resonant energy. 

Moreover, in the high voltage limit, eV ^> Erh, the 
differential Fano factor saturates at a constant value 2, 
i.e. independent on the relation between the conduc- 
tances g\ and <?2- This is a manifestation of the uncor- 
related charge transfer at these energies, as discussed in 
connection to Eq. (ful). Also in the limit <?2 ^> Q\, the 
Fano factor approaches 2, independent on voltage. This 
is again a manifestation of the uncorrelated pair transfer, 
and can be seen directly from the generating functions in 
Eq. (|l5|). From this we conclude that the differential 
Fano factor is thus the relevant quantity for studying the 
energy dependence of the charge transfer mechanism. 



IV. CONCLUSIONS 

In conclusion, we have investigated the energy 
dependent full counting statistics of a diffusive 
normal- insulating- normal-insulating-superconducting 
(NIxN' I 2 S) junction. We have used the recently devel- 
oped circuit theory of full counting statistics, allowing us 
to access the full temperature and voltage dependence 
of the statistics. In general, the charge is transported 
via correlated transfer of pairs of electrons. However, 
in the case of strongly asymmetric contacts or energies 
much larger than the Thouless energy, the pair transfer 
is uncorrelated. The second cumulant, the noise, was 
studied in detail. It was found to depend strongly on 
voltage and temperature. For low temperatures and a 
junction resistance dominated by the tunnel barrier to 
the normal reservoir (ii), the noise shows a double-peak 
behavior at voltages eV ~ Exh, a signature of an 
ensemble averaged electron-hole resonance. 
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